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Modeling and New Equipment Definition for the 
Vibration Isolation Box Equipment System 


R.L. Sani 

Department of Chemical Engineering 
University of Colorado 
Boulder, CO 80309-0424 


Our MSAD-funded research project is to provide numerical modeling support for the VIBES 
(Vibration Isolation Box Experiment System) which is an IML2 flight experiment being built 
by the Japanese research team of Dr. H. Azuma of the Japanese National Aerospace 
Laboratory. During this reporting period, the following have been accomplished: 

1. A Semi-Consistent Mass Finite Element Projection Algorithm For 2D And 3D 
Boussinesq Flows Has Been Implemented On Sun, HP And Cray Platforms. 

a. The algorithm has better phase speed accuracy than similar finite difference or 
lumped mass finite element algorithms, an attribute which is essential for 
addressing realistic g-jitter effects as well as convectively-dominated transient 
systems. 

2. The Projection Algorithm Has Been Benchmarked Against Solutions Generated Via 
The Commercial Code FIDAP. 

a. The algorithm appears to be accurate as well as computationally efficient. 

3. Optimization And Potential Parallelization Studies Are Underway. 

a. Our implementation to date has focused on execution of the basic algorithm with 
at most a concern for vectorization. 

4. The Initial Time-varying Gravity Boussinesq Flow Simulation Is Being Set-up. 

a. The mesh is being designed and the input file is being generated. Some 
preliminary "small mesh" cases will be attempted on our HP9000/735 while our 
request to MSAD for supercomputing resources is being addressed. 

5. The Japanese Research Team For VIBES Was Visited, The Current Set-up And Status 
Of The Physical Experiment Obtained And Ongoing E-Mail Communication Link 
Established. 



REQUEST FOR NASA SCIENTIFIC 
COMPUTING RESOURCES 


Robert L. Sani 
University of Colorado 
Department of Chemical Engineering 
Boulder, Colorado 80309-0424 


Our MSAD-funded research is to provide numerical modeling support for the 
VIBES (Vibration Isolation Box Experiment System) which is an IML2 flight experiment 
being built by the Japanese research team of Dr. H. Azuma of the Japanese National 
Aerospace Laboratory. This request of scientific computing resources is made to support 
flow and transport modeling for ground base pre-flight experimental design and 
validation experiments, post-flight data analysis as well as an assessment of the potential 
of using data and hardware for future flight experiments. In order to accomplish these 
goals, the transient Navier-Stokes and/or Boussinesq equations in three space dimensions 
including a complicated time varying gravitational body force must be solved. A new 
semi-consistent mass projection algorithm has been developed to generate time-accurate 
solutions. This algorithm has better phase speed accuracy than similar finite difference or 
lumped mass finite element algorithms and appears to be very efficient computationally, 
both necessary attributes for addressing realistic g-jitter effects as well as transient 
systems dominated by convective transport. (A brief description of the project and 
algorithm are in Appendix 1.) The development and preliminary benchmarking of the 
algorithm has been done on our HP9000 workstation, but access to a supercomputer 
platform is essential to support the very large 3D transient simulations especially 
necessary in the post-flight data assessment and in new experiment assessment; such 
simulations even when theoretically possible on a workstation are practically intractable 
because of the amount of computational time required. The algorithm has been designed 
such that the porting to a supercomputer such as a CRAY Y-MP will be very easy as will 
be the efficient use of its vectorization capability. Parallelization is an issue which we are 
currently addressing; the possibility of utilizing the efficient NASA parallel-vector direct 
solution algorithm due to Storaasli et. al. is already incorporated. Other parallelization 
via iterative solvers will eventually be addressed but only aggressively at a later date 
because of the time-line constraint of our MSAD project. The algorithm on the CRAY 
Y-MP is dynamically dimensioned and for the types of 3D meshes which seem to be 
appropriate for production runs, we will require approximately 64 megawords and an 
estimated 100 hours of computing time for algorithm porting enhancement, optimization 
and the production runs. The latter is difficult to estimate because of our desire to use an 
actual filtered g-jitter signature in the post-flight stage of the project. 



APPENDIX 1 



PROJECT DESCRIPTION 


A. Introduction 

In a terrestrial environment, natural convection and concomitantly the heat and mass 
transport processes intimately associated with it can sometimes have a deleterious effect on a 
process; for example, in materials processing where it, in conjunction with other effects such as 
walls or moving interfaces, can lead to radial segregation and dopant striations in the product 
material. In view of this, there has been a desire to use the space environment in order to reduce 
natural convection and other gravitational effects as well as make containerless processing more 
viable. In many cases, the ensuing dominance of surface -tension-driven flows in a micro-gravity 
environment has led to a significant number of similar as well as new problems. To further 
complicate matters, the micro-gravity environment available is not quiescent but is subjected to 
significant background vibrations generated by aerodynamic, structural and machinery vibrations, 
crew motion, etc. Such g-jitter can be relatively random in orientation and attain significant 
magnitudes such as 10- 2 g o (g 0 = earth g) experienced during thruster firing on the D1 mission to 
the less significant 10- 5 g o experienced on the SPAR-X free flyer. There is a growing list of 
observations and data analyses that demonstrate the existence of significant g-jitter episodes and 
the potential for having very deleterious effects on many proposed flight experiments as well as 
negate the potential of a micro-gravity environment. In addition, numerical and analytical 
modeling studies suggest that g-jitter effects could seriously affect many micro-gravity 
experiments as well as possibly space commercialization ventures, especially in the materials 
science area. An interesting potential solution to this problem in the micro-g environment 
available on shuttle flights or planned on Space Station Freedom is the use of vibration isolation 
for the experiments and processes which require it. The assessment of such an apparatus is one 
of the main thrusts of the research proposed herein. 

The experiment VIBES (Vibration Isolation Box Experiment System) is an IML2 flight 
experiment being designed by a Japanese research team under the direction of Dr. Hisao Azuma 
of the Japanese National Aerospace Laboratory. Its basic objective is to assess the performance 
of a vibration-isolation device in conjunction with typical micro-g fluids experiments. The 
device has been specifically designed for the low frequency range which appears to be most 
deleterious to many micro-g experiments, especially those material science experiments dealing 
with coupled flow and transport processes and could provide a relatively simple and economical 
means of dealing with such g-jitter problems. The IML2 flight experiment will contain two 
experimental units: 

1. Convection Diffusion Unit (CDU) 

2. Thermal-Driven Flow Unit (TDFU) 

The CDU experiment is the one of interest herein; the main objectives of this experiment 
are to observe natural convection and diffusive transport in a micro-g environment and to 
observe the effect of a g-jitter with and without the vibration isolation due to the vibration 
isolation box. The semi-consistent mass projection method algorithm summarized hereafter has 
been inplemented to provide flow and transport modeling capability for ground-based experiment 



design and validation and post-flight data analysis as well as an assessment of the potential of 
using the data and hardware for a future flight experiment. It is designed to solve the transient 
Navier-Stokes, or Boussinesq equations in two or three space dimensions. Time accurate 
solutions are obtained by using a semi-consistent mass projection algorithm [Gresho (1990)]. 
Theoretically this algorithm has better phase speed accuracy than similar finite difference or 
lumped mass finite element algorithms, a property which is very desirable for the g-jitter 
simulations being addressed as well as other convective transport dominated systems. 

CONTINUUM MODEL 

In terms of the Boussinesq equations, the appropriate continuum equations for the 
velocity (u ), temperature (T), and kinematic pressure ( P ) are: 

— + uVu = -V/ ) +vV : u + f , in 12 (la) 

dr 

V u = 0 , in Q (lb) 


pC p (^ + u VT) = kV 2 T+<? , in Q 

9 t 

The vector / represents body forces, for example, buoyancy. 
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f (P„-P(T)) 

buoyancy p e 

where g is the gravitational vector which in our case will be time-dependent, v is the kinematic 

viscosity, p (pj is the density (reference density), C p is the specific heat and k is the thermal 

conductivity. If the boundary of Q. is designated by 9 Q = T, © H , then typical appropriate 
boundary' conditions are specified velocity and temperature (Dirichlet): 

u = w , T = T b (Id) 

on f, and 'pseudo traction’ and specified flux conditions (Neumann) 
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on r 2 . Here n represents the outward normal direction, u n = u n is the normal component of 
velocity, x represents the corresponding tangential direction, w t = u x is the tangential 

component of velocity and F n and F t are the normal and tangential components of the specified 
boundary traction. (All combinations of BC's (Equations lc and Id) are possible in three- 
dimensional simulations)) 

To complete the problem specification, initial conditions must be specified: 

u(x,0) = u 0 (x) (If) 

T(x,0) = T 0 (x) 

The initial velocity field is required to satisfy the following conditions: 

n u 0 = n-w(x.O) on T, and (lg) 

V u 0 =0 inQ (lh) 

in order that a solution exists [Gresho and Sani (1987)]. Should T, = global mass 

conservation 


|n w(x, f) = 0 (li) 

an 


is an additional solvability constraint. 

The continuum version of the projection technique is a finite element discretization (with some 
modifications) of those proposed by Van Kan (1986) and Bell and Marcus (1990) as described in 
Gresho (1990). The algorithm is: 

0) Given u 0 with V • u 0 = 0 and P 0 , 

1) Solve for u and T, with u 0 = u 0 and T = T 0 at r= 0 from 


— + u • Vu - vV 2 u = f - V/l , in Q 

3 1 
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u = w and T = T„ , on T, , 
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for 0 < t < t v where [0,/J is the time span between two projections. Since the intermediate 

velocity u is not divergence free, solve the following least squares problem at projection time T : 
Extremize: 


G(v,cp)sji(v-u) 7 ‘(v-u)-<p r Vv (3a) 

Subject to appropriate boundary conditions; 

i.e., set the first variation equal to zero to obtain (after an integration by parts): 



(v-u) + V<p = 0 and 

(3b) 


V - v = 0 in Q , 

(3c) 

with 

n v = n u on r, and (p = ) + /J,] on T 2 


or in 

practice. 


2) 

Solve for (p from 



V : cp = V u in Q , 

(3d) 

subject to 



— = 0 on T, and 

an 

(3e) 


= onr,. 

(30 

3) 

compute v = u( T)- V<p in Q 

(4) 

4) 

Estimate a new pressure: P(T)= P 0 + 2(p / r, in Q 

(5) 



5) Report v and P; then set r=0, P 0 = P(t ] ), u 0 = v in £2 and go to step 1 . 


DISCRETIZATION 

The continuum formulation is discretized in space via a Galerkin finite element method. For 
efficiency, the finite element representation is restricted to four node (in 2D) or 8 node (in 3D) 
elements for velocity and temperature and a piecewise constant pressure approximation. Even 
though these element combinations are degenerate, i.e., the velocity-pressure basis function 
combination possesses a null space, it is easily accommodated [Sani, Gresho, Lee and Griffiths 
(1981)]. The use of a basic Galerkin finite element discretization technique automatically 

accommodates complex spatial domains and the complex boundary conditions on T, which are 
natural boundary conditions in the weak Galerkin finite element formulation. 

The continuous-in-time Galerkin finite element spatial discretization is then temporally 
discretized using the Leismann scheme. [Leismann and Frind (1989)] which is second-order 
accurate and unconditionally stable in the form implemented but a bit dissipative or the Gresho 
and Chan (1990) scheme which is similar but less dissipative. The discrete system for each time 
step of the Leismann scheme is: 

<— M + D+iD TO ) i u,-(-!-M-V-iD ra) ) 1 u a -MMi(C 1 P,-t(T)) . i-1. 2, 3 <«*) 

At 2 At 2 


( ±M + D + i-D BTD )T = (-J-M-V-^D BTD )T 0 +q 0 , 
At 2 At 2 

and 


(C T M['C)cp = C T u 

for tp needed to project the intermediate velocity u to the weakly divergence free subspace, i.e., 
the final, weakly divergence free velocity u after a time step is 

u = u - Mj^ 1 Cep (8) 

and the corresponding pressure field is: 

P = P 0 + 2tp / At (9) 


The various element matrices M, D, D btd V and M L are defined in Gresho and Chan 
(1990). The two most noteworthy features of this discretization are: 



1. The consistent mass matrix M is used in the conservation of momentum and energy 
equations (6a, b) but lumped mass matrix M L in the solenoidal velocity constraint equation 
(7). 

2. The discretized system of uncoupled algebraic equations (6a, i.e. 1,2,3; 6b) are linear 
symmetric positive definite systems. 

Feature (a) leads to phase speed accuracy while maintaining computational efficiency. Feature 
(b) allows the utilization of computationally efficient interative solvers such as pre- conditioned 
conjugate gradient solvers or, for example, the PVS-solver [Storaasli et. al.] (1990) which is a 
vectorized-parallelized efficient direct solution technique. 

IMPLEMENTATION 

The algorithm has been implemented on HP9000/735 and SUN/ 10 workstations and has 
been benchmarked against solutions to both Navier-Stokes and Boussinesq flows generated via 
the commercial code FIDAP [Engelman (1993)]. The current version utilizes primarily a pre- 
conditional conjugate gradient solver and a constant time step scheme with a variable time step 
version under development. The algorithm has been structured for easy vectorization and the 
possibility of also addressing the issue of parallelization, via the Storaasli algorithm initially and 
other iterative techniques eventually. 
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